Energy versus Angular Momentum in Black Hole Binaries 
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Using accurate numerical relativity simulations of (nonspinning) black-hole binaries with mass 
ratios 1 : 1, 2 : 1 and 3 : 1 we compute the gauge invariant relation between the (reduced) binding 
energy E and the (reduced) angular momentum j of the system. We show that the relation E(j) is 
an accurate diagnostic of the dynamics of a black-hole binary in a highly relativistic regime. By com- 
paring the numerical-relativity -E NR (j) curve with the predictions of several analytic approximation 
schemes, we find that, while the usual, non-resummed post-Newtonian-expanded _E PN (jf) relation 
exhibits large and growing deviations from _E NR (j), the prediction of the effective one-body for- 
malism, based purely on known analytical results (without any calibration to numerical relativity), 
agrees strikingly well with the numerical-relativity results. 

PACS numbers: 04.30.Db, 04.25.Nx, 95.30.Sf, 97.60.Lf 



Introduction. - A ground-based network of interfer- 
ometric gravitational wave (GW) detectors is currently 
being upgraded and is expected, thanks to an improved 
sensitivity, to detect, within a few years, the GW sig- 
nals emitted during the inspiral and merger of compact 
binaries. The realization of this exciting observational 
prospect depends, however, on our theoretical ability to 
accurately compute, within Einstein's theory of general 
relativity, the motion of compact binaries and its associ- 
ated GW emission. Recent developments have made it 
clear that the most efficient way to theoretically under- 
stand the late stages of the dynamics of compact bina- 
ries is to combine the knowledge coming from analytical 
relativity techniques -such as traditional post-Newtonian 
(PN) expansions (lnjl, or the newer effective-one-body 
(EOB) formalism |5|— lS|] , with the knowledge coming from 
numerical relativity (NR) simulations (see Q for a recent 
review). Here, we shall restrict our attention to binaries 
composed of two nonspinning black holes of masses mi 
and 77i2. Our technique can, however, be applied to more 
general systems. 

The aim of this Letter is to present how NR data can 
be used to explore, in a quite direct manner, the dy- 
namics of black-hole binaries, by computing the relation 
between the total energy, £ , of the binary system, and 
its total angular momentum, J. We compare the (gauge- 
invariant) relation £ (J) extracted from NR simulations 
to the corresponding analytical predictions from PN the- 
ory [l(|, and from EOB theory Q. We show that, dur- 
ing the inspiral (at least up to the last stable orbit) the 
gauge- invariant relation £{S) is essentially independent 
of the current uncertainties in the analytic modelling of 
the emitted gravitational waveform, and can therefore in- 
form us rather directly on the conservative dynamics of 
a black-hole binary. [This aspect of our work is akin to 
a recent study of periastron advance in black-hole bina- 
ries [ill.] 

Numerical relativity. - Our results are based on new, 
accurate numerical simulations of (nonspinning) black- 



hole binaries, which combine a 3 + 1 Cauchy-evolved 
spacetime (using a variant of the "BSSNOK" evolution 
system, with moving punctures and an extended wave 
zone [III, [l3|) with a Cauchy-characteristic extraction 
(CCE) technique [3 El. The initial data for the 3 + 1 
evolution are conformally flat, Bowen-York Cauchy data, 
with the initial position and linear momenta of the punc- 
tures determined from a 3PN-accuratc dynamical evolu- 
tion starting from a large initial separation (l6j . These 
initial data lead to orbits having an eccentricity e ~ 10~ 4 . 
The CCE technique yields unambiguous estimates of the 
waveforms at infinity, without the need to extrapolate 
data extracted at finite radii. Here, we consider three 
simulations with mass ratios q = milm\ equal to 1, 2 
and 3. The corresponding initial Arnowitt-Dcscr-Misncr 
(ADM) total energy, £q = £adm, total angular momen- 
tum, j7o = J7adm (oriented along the z axis), and eccen- 
tricity are given in Table I. 

We use these numerical simulations to compute the 
sequence of instantaneous values (at the retarded time t) , 
£{t), J{t), of the system energy and angular momentum 
during the inspiral, by using the laws of conservation of 
£ and J between the binary system and the emitted 
radiation. Namely, we compute 



nt) = ft-«t), 
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where the radiated energy and angular momentum, be- 
tween the initial (retarded) time to and time t, are com- 
puted from the multipolc moments of the NR (com- 
plex) "news function" at infinity (we generally use units 
such that G = c = 1): 
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TABLE I. Properties of the initial state of the NR simulations. 
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0.25 


1.5 x 10~ 4 


0.9905197 


0.9932560 


2 


2/9 


1.2 x 10~ 4 


0.9908980 


0.8559960 


3 


0.1875 


7.6 x 10~ 4 


0.9933905 


0.7675068 



E 



(4) 

Here is the NR multipolar metric waveform, 

Ni m (t) = dhi m {t)/dt and £ max = 8. We do not write here 
the expressions for the radiative losses of the other com- 
ponents J x , J y of J . We took them into account, though 
they turn out to have a negligible effect on the computa- 
tion of J NR (<). While AS™ only depends on the news 
function N(t) (which is a direct output of the CCE code), 
the angular momentum loss also depends on the metric 
waveform h(t). We computed (for each multipole) h(t) 
from ^>4(t) — dN/dt — d 2 h/dt 2 by the frequency-domain 
integration procedure of [171 ] (with a low-frequency cut- 
off ljq = 0.032/(mi + m2)). In contrast to most studies of 
gravitational waveforms, we consider here the full time 
development of the GW emission from the start of the 
NR simulation, i.e., we crucially take into account the 
losses associated with the "junk radiation" , viz the ini- 
tial burst of radiation associated to the relaxation of the 
unphysical Bowcn- York-type initial data, before the radi- 
ation settles down to a quasi-stationary inspiral pattern. 

Finally, we replace the two i-parametrized series 
£ NR (t), j7 NR (i) by the corresponding unparametrized 
curve £ NR (J). One example (for the mass ratio q = 1) 
of our computations of the relation £ (J7) is shown in 
Fig. 1. Here and below we work with the binding en- 
ergy per reduced mass, E = {£ — M)/fi, and the dimen- 
sionless rescaled angular momentum j = J/Mfi, where 
M = mi + ?ri2, t 1 = fn\rail (ra\ + 7712). Fig. 1 com- 
pares the NR relation P NR (j) to the predictions made 
by two different analytical formalisms: PN theory and 
EOB theory (as explained in detail below). The inset 
shows the very significant effect of the energy loss due to 
the junk radiation emitted at the beginning of the simu- 
lation. Note that j decreases during the inspiral. 

Post-Newtonian theory. - The gauge-invariant relation 
E(j) has been computed in PN theor y at the second post- 
Newtonian (2PN) approximation in [18[ and at the third 
post-Newtonian (3PN) one in (lb] . It has the structure 
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where c n (v) are polynomials (of order n) in the sym- 
metric mass ratio v = n/M = m\rri2l(rn\ + m^) 2 . [Sec 
Eq. (5.1) of [H], completed by for the final determina- 
tion of the 3PN dynamics, viz w s t a tic = 0.] The "Taylor" 
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FIG. 1. Equal-mass case: comparison between four E(j) 
curves. The standard "Taylor" PN curve shows the largest 
deviation from NR results, especially at low j's, while the 
two (adiabatic and nonadiabatic) 3PN-accurate, non-NR- 
calibrated EOB curves agree remarkably well with the NR 
one. 



(i.e., nonresummed) E PN (j) curve is shown in Fig. [T] (for 
q = 1) as a dashed line. 

Effective- one-body theory. - The EOB formalism 
maps the conservative dynamics of a two-body system 
onto the dynamics of one body of mass fi in a sta- 
tionary and spherically symmetric "effective" metric, 
ds 2 



-A(r; v)dt 2 + (A(r; u)D(r; v^dr 2 + r 2 {d6 2 + 
sva. 2 9dip 2 ). The EOB potentials A and D have been 
computed at the 2PN approximation in and at 
the 3PN approximation in Q (at 3PN one must com- 
plete the geodesic dynamics by terms, Q(p), quartic 
in momenta). Here, we use the 3PN-accurate ver- 
sion of the EOB Hamiltonian, as defined in 2000 0] 
(with astatic = Q), i.e., with the effective- metric 
potentials D(u) = 1 + dvu 2 + (52 — Qv)vu z , and 
A(u) ee P3 1 [l - 2u + 2vu 3 + (f - §±7r 2 ) vu^\ , where 
u ee GM/(c 2 r), and where P3 1 denotes constructing a 
(1, 3) Padc approximant, so that A(u) is a rational func- 
tion of u of the form (1 + ri\u)/{l + d\U + d^u 2 + d 3 u 3 ). 
In addition to the Hamiltonian dynamics defined by 
A(u), D(u) (and Q(u,p)), the EOB formalism defines 
a radiation-reaction force J- v . Here, we use the "newly 
resummed" radiation reaction defined by (l9l [20l | , with 
3 +2 -PN accurate Taylor pf m 's, and without incorporat- 
ing any "next-to-quasi-circular" (NQC) correction factor. 
The main point is that the resulting radiation-reaction- 
driven EOB dynamics uses only information that has 
long been analytically known, and does not rely on 
any information deduced from comparing EOB wave- 
forms to NR waveforms. The resulting (nonadiabatic) 
3PN-accurate, radiation-reaction driven EOB dynamics 
leads to the curve E EOB3rN (j) shown in Fig. [T]as a red 
solid line. In addition, we also show the adiabatic EOB 
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FIG. 2. Differences between seven E x (j) curves and 
E EOB3PN (j), for the three mass ratios considered. From top 
to bottom the labelling is: X = PN, EOB"pN Q °, EOB^? N c , 
NR, EOBspn (baseline), EOB^ N c and EOB3PN batic . While 
the PN curve exhibits the largest deviations, all EOB curves 
remain close to the NR one during the full inspiral, especially 
the 3PN-accurate, non-NR- calibrated one. 



i?EOB 3PN (J) curve defined by considering the sequence 
of minima in r (for a fixed j) of the (3PN- accurate) 
EOB Hamiltonian F EOB3PN (r, j). This adiabatic curve 
only depends on the potential A(u) and has a cusp at 
the last stable orbit (LSO), j'lso- The vertical distance 
£Eob 3PN (j) _ £ EOBS bstlc ^ essentially represents the 
kinetic energy linked to the (slow) inspiralling radial mo- 
tion. 

Results of the triple comparison NR-PN-EOB. - Fig. Q] 
already exhibits several of the new results of our study: 
(i) The NR E(j) curve starts at large j's (i.e., large radial 
separations) close to the PN- predicted E PN (j) curve, but 
then visibly deviates more and more from it during the 
inspiral [conventionally ending at the adiabatic (-EOB- 
defined) LSO, marked by a dashed vertical line], (ii) By 
contrast the NR E(j) curve is so close, on the scale of 
Fig. [TJ to the (3PN-accurate, nonadiabatic) EOB predic- 
tion that their difference is barely visible not only during 
the inspiral, but also during the subsequent plunge. [The 
leftmost red vertical line in Fig. [TJ denotes the EOB "light 
ring" , viz the end of the analytical inspiral-plus-plungc 
dynamics, and the beginning of the EOB description of 
the merger and ringdown.] (iii) On the scale of Fig. [TJ one 
cannot see, during the inspiral, the difference between the 
two EOB curves (nonadiabatic versus adiabatic). (iv) In 
addition, when zooming on the beginning of the £ ,NR (j) 
curve (see inset), we find that, although it coincidcntally 



starts near the PN curve, it emits exactly the amount of 
junk radiation required to relax to the EOB prediction. 
When considering the mass ratios q = 2 and q = 3, we 
obtained close analogs of Fig. [TJ which exhibit exactly 
the same results (i)— (iv). 

In order to refine and quantify these results, we hence- 
forth zoom on the small deviations between the var- 
ious E(j) curves by using as horizontal baseline the 
(nonadiabatic, 3PN-accurate) EOB curve, i.e., by plot- 
ting the differences E x (j) - £ EOBaPN (j), where the label 
X denotes either NR, PN, EOB^™ batic , or other EOB 
curves defined below. When focussing on the inspiral 
dynamics (above the LSO), this leads to NR-EOB dif- 
ferences of order 10 -4 , i.e., 300 times smaller than the 
~ 3 x 10 -2 change in the absolute value of E during the 
inspiral, and 50 times smaller than the PN-NR differ- 
ence ~ 5 x 10~ 3 at the LSO. To discuss the meaning of 
the small NR-EOB differences, it is important to esti- 
mate the error attached to the NR E NR (j) curve. We 
estimate an error on E m (j) by measuring the effect 
of changing, in turn, all the NR elements entering the 
computation of the losses Eqs. (0-(j4|): (i) we replaced 
the CCE news by either the time integral of the cur- 
vature waveform ^(t) = dN{t)/dt extracted at a large 
radius in the 3 + 1 code, or a Regge-Wheeler-Zerilli met- 
ric waveform output by the latter code; (ii) we reduced 
the maximum multipolar order £ max used in the sums in 
Eqs. ©-dU) from the default value £ max = 8 to £ max = 7 
and £ roax = 6; (iii) we varied the low-frequency cut- 
off Mujq used in the frequency-domain computation of 
h hn (t) from ^ e 4 m (t) [13 between about 0.01 and 0.04; 

(iv) we computed hg m (t) from N^ m (t) instead of ^^"(i); 

(v) we explored the sensitivity to changes of the initial 
integration time t in Eqs. J3])-(J4|); (vi) we replaced the 
high resolution NR data used as a baseline by medium 
resolution ones. 

Adding the effect of all these changes, and focussing 
on the crucial change in the energy loss A£j un k linked 
to the initial burst of junk radiation, leads to a con- 
servative error bar around £ ,NR (j) indicated by a gray- 
shaded region in Fig. [3J In that figure, we plot the dif- 
ferences E x (j) - E EOB3 ™(j) for q = 1, 2 or 3, and for 
six different labels X: NR (solid, thick, black curve), PN 
(upper, dashed blue curve), EOB3p^ batIC (lower, dash- 
dotted magenta curve), EOB^j (black, dashed curve, 
just below the baseline) , EOB^j (upper solid red curve) 
and EOB™n QC (dashed red curve, close to the previous 
one). Here, as above, the EOB baseline i<; EOB 3PN (corre- 
sponding to the horizontal axis), as well as its adiabatic, 
EOB a P ^ batic , and NQC-completed, EOB^ , avatars, 
use the 3PN-accurate EOB potentials of Q. [EOB^ 
is defined according to the methods introduced in |2Cf by 
adding a factor /^^ c (01,02) in the £ = m = 2 mode, 
tuned to the maximum of the NR modulus.] Finally, 
£ EOB S° and £ eob ™p°n QC use the NR-calibrated, 5PN 
potential A^p-^(u) = [A (u) + va^u 5 + i>a§u & \, 
for (05,06) = (—6.3722,50) [which lies in the middle 



4 



of the "good region" of Ref. 0], either with (NQC) 
or without (wo NQC) NQC corrections. Figure [2] al- 
lows us to refine and strengthen the conclusions drawn 
above from Fig.[TJ namely: (i) The PN- predicted E PN (j) 
curve disagrees strongly with the NR results; (ii) The 
3PN-accurate nonadiabatic EOB curve, E EOB3FN (i.e., 
the horizontal baseline) is remarkably close to the NR 
results during the entire inspiral, with deviations that 
are smaller than the "2<t" level, (iii) The inclusion of 
nonadiabatic effects is important in continuing to ensure 
this agreement during the late inspiral (see the difference 
E EOB^ t,c - £EOB3pn). ( iv ) The i nc i us ion of the NR- 
fitted NQC correction has a negligible effect during the 
inspiral: £ woNQ c - ^nqc < 2 x 1CT 5 . (v) The EOB 
predictions based on the NR-calibrated, 5PN potential 
^5Pn(u) of Ref . [2(| (with or without NQC corrections), 
are less close (especially for q = 1 and 3) to the NR result 
than the purely analytical 3PN-accurate EOB prediction. 
We verified that the same conclusion holds for the NR- 
calibrated 5PN EOB potential suggested in (2lj . 

Summary. - We showed how to combine the knowledge 
of the initial (ADM) energy and angular momentum of 
a black-hole binary with accurate NR computations of 
its subsequent GW emission (including the initial burst 
of junk radiation), to derive the relation between the 
rescalcd binding energy E = (£ — M)//z and the rescalcd 
angular momentum j = J /(Mfx). Though the relation 
E(j) does include nonadiabatic effects (linked to the ra- 
dial kinetic energy during the inspiral, and thereby to the 
radiation reaction T^) we have verified that the analytic 
uncertainties in the description of were essentially 
negligible during the inspiral, down to, at least, the LSO. 



This makes the NR-acquired knowledge of the E(j) curve 
an accurate diagnostic of the conservative dynamics of a 
black-hole binary in a highly relativistic regime. By com- 
paring E (j) to various analytic descriptions of binary 
dynamics, we found that, while the usual, nonresummed 
3PN-expanded relation i? PN (j) exhibits large and grow- 
ing deviations with respect to i? NR (j), the EOB formal- 
ism, based purely on known analytical results (without 
NR calibration) predicts a relation E EOB (j) which is re- 
markably close to E NR (j). We found that the various ex- 
istent NR-calibrated EOB formalisms fare somewhat less 
well than the purely analytic 3PN-accurate EOB formal- 
ism in agreeing with the NR results. This clearly shows 
that the NR curve i? NR (j) contains valuable information 
about the conservative dynamics of the binary during the 
inspiral that can usefully complement the information 
contained in the waveform (which mixes in an intimate 
manner dynamical and radiative effects). We leave to fu- 
ture work a detailed discussion of the information that 
can be extracted from E NR (j). 
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